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■ ABSTRACT 

o : 

O ; 

CN | Observations of the Lya forest at z ~ 3 reveal an average metallicity Z ~ 10 2 Zq. 

The high-redshift supernovae that polluted the IGM also accelerated relativistic elec- 

h-^ ■ trons. Since the energy density of the CMB oc (1 + z) 4 , at high redshift these electrons 

cool via inverse Compton scattering. Thus, the first star clusters emit X-rays. Unlike 
stellar UV ionizing photons, these X-rays can escape easily from their host galaxies. 
. This has a number of important physical consequences: (i) Due to their large mean free 

path, these X-rays can quickly establish a universal ionizing background and partially 

\Q \ reionize the universe in a gradual, homogeneous fashion. If X-rays formed the domi- 

nant ionizing background, the universe would have more closely resembled a single-phase 
medium, rather than a two-phase medium, (ii) X-rays can reheat the universe to higher 
temperatures than possible with UV radiation, (iii) X-rays counter the tendency of UV 

■ radiation to photo-dissociate H_2, an important coolant in the early universe, by pro- 
moting gas phase H2 formation. The X-ray production efficiency is calibrated to local 
observations of starburst galaxies, which imply that ~ 10% of the supernova energy is 
converted to X-rays. While direct detection of sources in X-ray emission is difficult, 
the presence of relativistic electrons at high redshift and thus a minimal level of X-ray 
emission may be inferred by synchrotron emission observations with the Square Kilo- 
meter Array. These sources may constitute a significant fraction of the unresolved hard 
X-ray background, and can account for both the shape and amplitude of the gamma- 
ray background. This paper discusses the existence and observability of high-redshift 
X-ray sources, while a companion paper models the detailed reionization physics and 
chemistry. 



1. Introduction 

While the theoretical literature on the epoch of reionization is large and increasing rapidly, 
our empirical knowledge of this period in the history of the universe is scant and may be succintly 
summarized: (i) The universe is likely to have been reionized in the period 5.8 < z r < 35, where the 



1 Current address: Theoretical Astrophysics, Mail Code 130-33, California Institute of Technology, Pasadena, CA 
91125; peng@tapir.caltech.edu 



- 2 - 



lower bound arises from the lack of Gunn- Peterson absorption in the spectra of high redshift quasars 
(Fan et al 2000), and the upper bound comes from the observation of small scale power in Cosmic 
Microwave Background anistropies (Griffiths et al 1999). In CDM cosmologies, non-linear objects 
above the cosmological Jeans mass 1O 4_5 M first collapse during this period, (ii) The presence of 
metals in the Lya forest implies that significant star formation took place at high redshift (Songailia 
& Cowie 1996, Songailia 1997). (iii) COBE constraints on the Compton y-distortion of the CMB 
(Wright et al, 1994, Fixsen et al 1996) implies that the IGM was not heated to high temperatures. 
This means it is unlikely that the hydrogen and helium in the IGM were collisionally reionized. 

Thus, the current state of observations is consistent with scenarios in which the universe was 
reionized by an early generation of stars or quasars. While we do not know whether stars or quasars 
were the dominant source of ionizing photons, the observational and theoretical case for quasars is 
somewhat more uncertain. Extrapolation of empirical quasar luminosity functions to high redshift 
do not yield enough ionizing photons to maintain the observed lack of Gunn- Peterson absorption at 
z ~ 5 (Madau, Haardt & Rees 1999). The requisite steepening of the faint end slope at high redshift 
necessary to boost ionizing photon production is constrained by the lack of red, pointdike sources 
in the Hubble Deep Field (Haiman, Madau & Loeb 1999); the authors find that AGN formation 
must have been suppressed in halos with v c < 50 — 75kms _1 . We do not have a sufficiently firm 
understanding of the formation and fueling of supermassive black holes to assert on theoretical 
grounds that AGNs must have been present at high redshift. On the other hand, a minimal level 
of high-redshift star formation is guaranteed by the observed metal pollution of the IGM. 

Theoretical scenarios in which stars or quasars figure predominantly have been calculated in 
detail. Our ignorance of the efficiency of gas fragmentation, and star/black hole formation as a 
function of halo mass, make the prediction of observable differences between these two scenarios very 
uncertain. Indeed, if one normalizes assumed emissivities to a fixed reionization epoch, differences 
between the two scenarios boil down to: (i) stars result in supernovae, which inject dust, metals, 
and entropy into the host galaxy and surrounding IGM, which affects subsequent chemistry and 
cooling, (ii) quasars have a significantly harder spectrum than stars. In particular, they produce 
X-rays. 

In this paper, I emphasize a hitherto neglected fact: high redshift supernova also produce X- 
rays, both by thermal emission from the hot supernova remnant, and inverse Compton scattering of 
soft photons by relativistic electrons accelerated by the supernova. Considerable X-ray emission is 
already observed in starburst galaxies at low redshift (e.g., Rephaeli et al 1995), and the efficiency 
of most proposed X-ray production mechanisms should increase with redshift (e.g., explosions 
take place in a denser medium at high redshift, hardening expected thermal emission; inverse 
Compton scattering becomes more efficient since the CMB provides a ready supply of soft photons 
Ucmb oc (1 + z) A ). Thus, the SED of high-redshift star forming regions is considerably harder 
than has been previously assumed. This blurs the distinction between stellar /quasar reionization 
scenarios, and has a number of important physical consequences: 
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• Escape fraction The escape fraction of UV ionizing photons in the local universe is small, 
~ 3 — 6% (Leitherer et al 1995, Bland-Hawthorne & Maloney 1999, Dove, Shull & Ferrara 
2000) and is expected to decrease with redshift (Wood & Loeb 1999, Ricotti & Shull 1999). 
On the other hand, X-rays can escape freely from the host galaxy. Thus, processing by the 
host ISM may imply that the universe was reionized by a significantly harder spectrum than 
previously assumed. 

• Reionization topology Photons from stellar spectra have a short mean free path and thus a 
sharply defined ionization front. This fact gives rise to the conventional picture of expanding 
HII bubbles embedded in the neutral IGM. The spectra of quasars is significantly harder and 
exert an influence over a larger distance (which is why it is much more difficult to perform 
numerical simulations of reionization by quasars; see Gnedin 1999). Nonetheless, for quasars, 
uL v oc z^~ ' 8 (Zheng et al 1997, although note that the observations were only in the radio- 
quiet AGN subsample at energies up to 2.6 Ry) and most of the energy for ionization lies just 
above the Lyman edge. Thus, there is still a sharply defined HII region and a thin ionization 
front where the ionization fraction drops sharply. By contrast, for the inverse Compton case 
vL v ~ const, there is equal power per logarithmic interval, and thus there is no preferred 
energy scale. In particular, there is no preferred scale for the mean free path of ionizing 
photons. When the universe is largely neutral, it is optically thick even to hard photons and 
all photons with energies E < E^^ = 1.5 (^^)°' 5 x^keV (where xhi is the mean neutral 
fraction) are absorbed across a Hubble volume. While the (more numerous) soft photons 
can only travel a short distance before ionizing neutral HI and Hel, the (less numerous but 
more energetic) hard photons will be able to travel further and ionize an equivalent number of 
photons by secondary ionizations. Thus, even if sources are distributed very inhomogeneously, 
reionization will be a fairly homogeneous event, with a largely uniform ionizing background 
and fluctuations in ionization fraction determined mainly by gas clumping. Instead of an 
two-phase medium in whose HII filling fraction increases with time, the early IGM may have 
been a single phase medium whose ionization fraction increases with time. 

• Increased reheating A hard spectrum can reheat the IGM to considerably higher temper- 
atures than soft stellar spectra, both through photoionization heating and Compton heating. 
A soft spectrum loses thermal contact with the IGM once HI and Hel are completely ionized 
(at the mean IGM density, HI has a recombination time longer than the Hubble time for 
z < 10), and the gas cools adiabatically due to the expansion of the universe (Hui & Gnedin 
1997). By contrast, a hard spectrum can continually transfer large amounts of energy from 
the radiation field to the IGM by ionizing Hell, which recombines rapidly (Miralda-Escude & 
Rees 1994). This feedback mechanism is important in increasing the Jeans mass, a proposed 
mechanism for preventing excessive cooling and star formation at high redshift (e.g., Prunet 
& Blanchard 1999). The higher IGM temperatures may also explain why observed Lya forest 
line widths are commonly in excess of that predicted by numerical simulations (Theuns et al 
1999, Ricotti et al 2000). 
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• Early universe chemistry H2 is an extremely important coolant in the metal-free early 
universe. While the neutral IGM is optically thick to UV ionizing photons, it is optically 
thin to photons longward of the Lyman limit (except at wavelengths corresponding to higher 
order hydrogen Lyman resonance lines, as well as H2 resonance lines). In particular, photons 
in the 11.2-13.6 eV range quickly establish a soft UV background which photodissociates H 2 
via the Solomon process, shutting down subsequent star formation (Haiman, Rees & Loeb 
1997, Cicardi, Ferrara & Abel 1998), unless the H2 opacity is sufficient to reduce the photo- 
dissociation rate (Ricotti, Gnedin & Shull 2000). If X-rays are present in the early universe, 
they can counter this H2 destruction. The IGM is also optically thin to X-rays, which can 
penetrate dense clouds of gas and promote gas phase H2 formation H + e~ — > H™ + 7 and 
H + H — > H2 + e~ by increasing the abundance of free electrons. Haiman, Abel & Rees 
(1999) show that if quasars were the dominant ionizing sources in the early universe, gas 
cooling and thus star formation can continue unabated. In Paper II, I show that in fact even 
if only stars were present, a self-consistent treatment of the stellar SED incorporating X-rays 
produced by supernovae favours H2 formation over destruction in dense regions. 

In this paper, I study the emission mechanisms and observational signatures of X-ray bright 
star clusters at high redshift. In Paper II (Oh 2000a), I address the changes in reionization topology, 
reheating and early universe chemistry mentioned above due to these X-rays. 

In all numerical estimates, I assume a background cosmology given by the 'concordance' values 
of Ostriker & Steinhardt (1995): {Q m , Sl A , Sl b , h, a 8h -i , n) = (0.35,0.65,0.04,0.65,0.87,0.96). This 
corresponds to Q b h 2 = 0.017, compared with n b h 2 = 0.020 ± 0.002 (95%c.l.) (Buries, Nollett & 
Turner 2000), and n b h 2 = 0.0205 ± 0.0018 (O'Meara et al 2000) from Big Bang Nucleosynthesis, 
and the significantly higher values 0.022 < Q b h 2 < 0.040(95%c.l.) (Tegmark & Zaldarriaga 2000) 
preferred by recent CMB anisotropy data, such as Boomerang and Maxima. 

2. Emission mechanisms 

2.1. Star formation at z > 3 

What fraction of present day stars formed at high redshift? Estimates of the comoving star 
formation rate as a function of redshift (Madau et al, 1996) should be regarded as lower bounds, 
particularly at high redshift, due to the unknown effects of dust extinction, and star formation 
in faint systems below the survey detection threshold. Indeed, Lyman break survey results for 
3.6 < z < 4.5 (Steidel et al 1999) suggest that after correction for dust extinction, the comoving star 
formation rate for z > 1 is constant, rather than falling sharply as previously believed. Furthermore, 
in recent years compelling evidence has emerged that the majority of stars in ellipticals and bulges 
formed at high redshift, z > 3. This comes from the tightness of correlations between various 
global properties of ellipticals which indicate a very small age dispersion and thus a high redshift of 
formation, unless their formation was synchronized to an implausible degree. The evidence includes 
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the tightness of the fundamental plane and color magnitude relations for ellipticals, and the modest 
shift in zero-points for these relations with redshift (Renzini 1998 and references therein). Since 
spheroids contain ~ 30% of all stars in the local universe (King & Ellis 1985, Schechter & Dressier 
1987), this would imply that ~ 30% of all stars have formed at z > 3. Since ~ 20% of baryons have 
been processed into stars by the present day (Fukugita, Hogan & Peebles 1998), this implies that 
~ 0.2 x 0.3 ~ 6% of baryons have been processed into stars by z ~ 3. Assuming widespread and 
uniform enrichment and 1 M of metals per 100 Mq of stars formed, this translates into an IGM 
metallicity of 6 x 10" " 4 ~ 3 x 1O~ 2 Z . At z ~ 3, the metallicity of damped Lya systems appears 
to be ~ O.O5Z (Pettini et al 1997), in reasonable agreement. The observed metallicity of the Lya 
forest at z ~ 3, to which most models of reionization have been normalized, is between 10~ 2 Z & 
and 1O~ 3 Z (Songailia & Cowie 1996, Songailia 1997), which would imply that only 0.2 - 2% of 
present day stars formed at z > 3. However, its metallicity may be more representative of low 
density regions, rather than the mean cosmological metallicity (Cen & Ostriker 1999). Note that 
normalization of high redshift star formation to Lyman a forest metallicities assume efficient metal 
ejection (which underestimates star formation if a significant fraction of metals are retained) and 
mixing (which overestimates star formation if Lya lines are preferentially observed in overdense 
regions which are sites of star formation). In this context, it is worth mentioning claims that 
Lya lines with 10 13 ' 5 cm~ 2 < Nhi < 10 14 5 cm -2 reveal lower metallicities by a factor of 10 than 
clouds with Nhi > 10 14 - 5 cm~ 2 (Lu et al 1999). Ellison et al (2000) find no break in the power 
law column density distribution for C IV down to log N(C IV) =11.7, and Schaye et al (2000) 
detect O VI down to thi ~ 10 _1 in underdense gas, so it appears that metal pollution was fairly 
widespread. I regard Z ~ 10~ 3 — 2.5 x 1O~ 2 Z at z ~ 3 as a fairly firm bracket on the range 
of possibilities. For inverse Compton radiation, this corresponds to an energy release per IGM 
baryon of e SN ~ 10 ( io^z e ) (cn") e V (where e is the efficiency of conversion of supernova energy 
to X-rays), which is comparable to the energy release in stellar UV radiation for the low escape 
fractions expected, e stellar ~ 10 ( io-^z Q ) (o1f) e ^' where f esc is the escape fraction of ionizing 
photons from the source. 



2.2. X-ray emission in local starbursts 

Most models of reionization use population synthesis codes to estimate the spectral energy 
distribution of starbursts. However, there are many processes associated with star formation that 
generate UV and X-rays, beside stellar radiation: massive X-ray binaries, thermal emission from 
supernova remnants and hot gas in galactic halos and winds, inverse Compton scattering of soft 
photons by relativistic electrons produced in supernovae. Indeed, X-ray emission appears to be 
ubiquitous among starbursts (e.g., Rephaeli, Gruber, & Persic 1995), and starburst galaxies may 
account for a significant portion of the XRB (Bookbinder et al 1980, Rephaeli et al 1991, Moran, 
Lehnert & Helfand 1999). The X-ray emission from these processes, which hardens the spectrum of 
starbursts and changes both the topology and chemistry of reionization, has to date been neglected 
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in studies of the z > 5 universe. 

The X-ray luminosity of starbursts correlates well with other star formation indicators; for 
example, David et al (1992) find a roughly linear relation between Lfir and Lx- As a very 
rough empirical calibration, the starburst galaxies M82 & NGC 3256 observed with ROSAT 
and ASCA (Moran & Lehnert 1997, Moran, Lehnert & Helfand, 1999) follow the relations (af- 
ter correction for absorption): L x ,o.2-iokcV = 8 x W~ 4 L IR , and L x ,5kcY = 1-2 x 10 4 Lr j5 ghz 
(note that 5 keV flux density is relatively unaffected by photoelectric absorption or soft ther- 
mal emission). The starburst model of Leitherer and Heckman (1995) yields ~ Lfir ~ 
1.5 x 10 10 (SFR/1 M0yr _1 )L ; as a cross-check, the empirical relation for radio emission is (Condon 
1992) L R = 1.4 x lO 28 (i//GHz)- a (SFR/M yr- 1 )ergs~ 1 Hz~ 1 , where a ~ 0.8. Together I obtain: 



One should not regard this as more than a rough order of magnitude estimate; a large scatter is 
expected in this relation. By way of comparison, Rephaeli et al (1995) obtain from the mean of 51 
starbursts observed with Einstein and HEAO, Lv,2-30keV ~ 8 x 10 _3 L/#, an order of magnitude 
greater; and David et al (1992) obtain from a sample of 71 normal and starburst galaxies a ratio 
lower by about an order of magnitude, largely due to the inclusion of normal galaxies (this is 
consistent with an inverse Compton origin for X-rays, since normal galaxies have much lower 
radiation field energy densities and would not be expected to show significant inverse Compton 
emission) . 

The typical observed X-ray spectrum is a power-law, L u oc z^~ ' 8 , consistent with a non-thermal 
origin. An obscured AGN is not likely to be the source of these X-rays, as several observations 
suggest that the X-ray emission is powered primarily by massive stars. In NGC 3256, observations 
by ISO fail to detect high excitation emission lines (Rigopoulou et al 1996). In M82, the optical 
spectrum is HH-like (Kennicutt 1992), discrete nuclear radio sources are spatially resolved (Muxlow 
et al 1994), its nuclear X-ray emission is extended (Bregman et al 1995), and the expected broad 
Ha emission is not detected (Moran & Lehnert 1997); the HEX continuum and Fe-K line emission 
of NGC 253 as observed by BeppoSAX is extended (Cappi et al 1999). Moran & Lehnart (1997) 
and Moran, Lehnart, & Helfand (1999) have modelled the X-ray emission of M82 and NGC 3256, 
and find inverse-Compton emission to be the most likely mechanism, rather than an obscured 
AGN or massive X-ray binaries. The ratio between the observed radio and X-ray fluxes (note 
that Lx/L syn oc Ujr/Ub, where Um is the energy density of the infra-red radiation field, and Ub 
is the energy density of the magnetic field) is consistent with an inverse-Compton origin for the 
X-rays. Furthermore, the X-ray and radio emission have the same spectral slope, as is expected if 
both types of emission are non-thermal, arising from the same population of electrons. The X-ray 
luminosity is energetically consistent with a star formation origin. Assuming a Salpeter IMF and 
that each supernova explosion yields 10 51 erg in kinetic energy yields an energy injection rate into 
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the ISM: 




(2) 
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where e is the fraction of energy injected into relativistic electrons; consistency with the observed 
value, (1), implies that e ~ 10%. Hereafter, I shall use the empirical equation (1) as a fiducial 
conversion between X-ray luminosity and star formation rate. 

Is an acceleration efficiency of e ~ 10% reasonable? The acceleration mechanism for relativistic 
electrons is poorly understood. While first-order Fermi acceleration in shocks is widely accepted 
as the acceleration mechanism for cosmic rays (Blandford &: Eichler 1987, Jones & Ellison 1991), 
the electron acceleration is thought to be more problematic, due to the smaller electron gyroradius 
(which leads to greater difficulties in bouncing an electron back and forth across a shock of finite 
thickness), and the difficulty of initially boosting the electron to relativistic speeds, where Fermi 
acceleration can operate (Levinson 1994). From measurements of cosmic rays energy density it is 
inferred that ~ 10% of the supernova kinetic energy, or ~ 10 50 erg per explosion, is liberated as 
cosmic rays (Volk, Klein, & Wielebinski 1989), but the division between electrons and protons at 
the source is not known. Since the measured ratio of cosmic ray protons to electrons is ~ 75 (e.g., 
Gaisser 1990), it might well be that relativistic electrons only constitute e ~ 10~ 3 of the supernova 
energy budget. Thus, the reader should be cautioned that e ~ 0.1, which corresponds roughly equal 
energy division between protons and electrons (e/p=l), may be an overly optimistic estimate of 
the energy injection into relativistic electrons. Theoretical models of shock acceleration, in which 
e/p is a free parameter, often set e/p~ 1 — 5% for consistency with cosmic-ray experiments (Ellison 
& Reynolds 1991, Ellison et al 2000). However, this is somewhat model-dependent: in models 
where electrons are injected directly from the thermal pool, ~ 5% of the energy in the shock 
must go to non-thermal electrons in order to match gamma-ray observations (Bykov et al 2000). 
Furthermore, note that the observed cosmic-ray e/p ratio could equally well be the result of different 
transport processes and energy loss mechanisms for electrons and protons. In particular, cosmic 
ray electrons are subject to loss processes which operate on much longer timescales for cosmic ray 
protons (inverse Compton, synchrotron losses, etc); the cosmic ray flux at earth for electrons could 
arise from a much smaller effective volume than that for protons. At the source, the energy division 
between protons and electrons could range between 1 and 100. Perhaps the most reliable means of 
inferring the proton/electron energy division is by direct observations of supernova remnants. In 
modelling the observed production of gamma-rays in the supernova remnants IC 443 and 7 Cygni 
observed by the EGRET instrument on the Compton Gamma Ray Observatory, Gaisser, Protheroe 
& Stanev (1998) find that a proton to electron ratio of 3-5 gives the best fit to the observed spectra, 
implying e ~ 0.02 — 0.03. Similarly, in modelling the 7-ray flux from 2EG J1857+0118 associated 
with supernova remnant W44, de Jager & Mastichiadis (1997) find e ~ 0.09. They speculate that 
electron injection by the pulsar may be responsible for the increased electron energy content. Given 
the large uncertainties, henceforth I shall simply use the empirical relation (1). 

Thus, barring non-standard IMFs, type II detonation energies or alternate sources of relativistic 
electrons, the empirical relation (1) implies an acceleration efficiency e ~ 0.1 which is plausible but 
certainly lies at the upper limit of theoretical expectations. Another possibility is that the X-ray 
emission cannot be wholly attributed to inverse Compton emission alone (this assumption rests on 
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the arguments of Moran & Lehnart (1997) and Moran, Lehnart, & Helfand (1999) with regards 
to the slopes and relative intensities of the observed non-thermal radio and X-ray emission). The 
X-ray emission may instead be due to X-ray binaries, thermal emission from supernova remnants 
or starburst driven superwinds (e.g., see Natarajan & Almaini 2000). It should be noted that 
only soft X-rays are relevant for reionization, since the universe is optically thin to photons with 
energies E > E t hick = 1-5 (txt)° 5 keV (where xhi is the mean neutral fraction of the IGM). 
Since equation (1) is calibrated with the soft bands observed by ROSAT, this implies that even 
if e <C 0.1 and the observed X-rays are not predominantly due to inverse Compton emission, 
the importance of X-rays for reionization (in particular, for changing the topology, for increased 
reheating, and increased H2 production) may still hold. However, observational predictions which 
focus specifically on the inverse Compton mechanism (e.g., the gamma-ray background (section 
(3.2)), and detecting synchrotron emission with the SKA (section (3.4)) will no longer be valid. 
Since the acceleration efficiency is the most uncertain parameter in this paper, wherever relevant I 
insert the scaling factor (g^-) into numerical estimates. 

How does the X-ray luminosity compare with stellar UV ionizing radiation? Assuming a 
Salpeter IMF with solar metallicity, the Bruzual & Chariot (1999) population synthesis code yields 
an energy output of Lj on = 3.2 x 10 42 (SFR/1Mq yr^ 1 ) ergs -1 in ionizing photons, which translates 
into Njon = 10 53 (SFR/1M© yr -1 ) photons s -1 . However, note that most of these ionizing photons 
are absorbed locally with the ISM of the star cluster; the escape fraction of ionizing photons into 
the IGM is expected to be small. Leitherer et al (1995) have observed four starburst galaxies 
with the Hopkins Ultraviolet Telescope (HUT). Their analyis suggests an escape fraction of only 
3%, based on a comparison between the observed Lyman continuum flux and theoretical spectral 
energy distributions. For our own Galaxy, Dove, Shull &; Ferrara (2000) find an escape fraction 
for ionizing photons of 6% and 3% (for coeval and Gaussian star formation histories respectively) 
from OB associations in the Milky Way disk. Bland-Hawthorn & Maloney (1999) find an escape 
fraction of 6% is necessary for consistency with the observed Ha emission from the Magellanic 
stream and high velocity clouds. On the other hand, a recent composite spectrum of 29 Lyman 
break galaxies (LBGs) with redshifts (z) = 3.40 ± 0.09 shows significant detection of Lyman con- 
tinuum flux (Steidel, Pettini & Adelberger 2000); for typical stellar synthesis models, the observed 
flux ratio L(1500)/L(900)=4.6 ± 1.0 implies little or no photoelectric absorption. The fraction of 
900 A photons which escape, f esc ~ 15 — 20%, is modulated almost entirely by dust absorption. 
Nonetheless, the authors themselves stress this result should be treated as preliminary; the result 
could be due to a large number of uncertainties or selection effects, among them the fact that 
these galaxies were selected from the bluest quartile of LBGs. On the theoretical side, radiative 
transfer calculations by Woods & Loeb (1999) find that the escape fraction at z ~ 10 is < 1% for 
stars; calculations by Ricotti Sz Shull (1999) find that the escape fraction decreases strongly with 
increasing redshift and halo mass; for a 10 9 M Q halo at z = 9, the escape fraction is ~ 10 -3 (note 
that in the Ricotti & Shull (1999) models, the escape fraction rises towards low masses, and can be 
considerable for the halos with M < 1O 7 M . Since such halos have T v i r < 10 4 K, their contribution 
to reionization depends on whether H2 formation and cooling can take place despite photodissocia- 



-9- 



tive processes (Haiman, Rees & Loeb 1997, Cicardi, Ferrara & Abel 1998, Ricotti, Gnedin & Shull 
2000)). For low escape fractions, the energy release in UV photons is roughly comparable to that 
in inverse-Compton X-rays : 

L uv = 3 x 10 40 ( ( SFR T ) ergs" 1 (3) 

Note that to first order the ratio of stellar UV to inverse Compton X-rays is not sensitive to 
uncertainties in the IMF, as the same massive stars with M > 20M Q that dominate the Lyman 
continuum of a stellar population also explode as supernovae (however, see section (2.4) for some 
caveats) . 

Will X-ray emission still be efficient at high redshift? The following changes are expected 
to take place at high redshift: (i) The ISM is initially free of dust and metals, although rapid 
enrichment could occur on fairly short timescales (t ~ 10 6 — 10 7 yr). (ii) The average ISM density 
is significantly higher, rihalo ~ nol87r 2 (l + z) 3 = 0.02( ^j) 3 cm -3 in a halo and n^ sc ~ X~ 3 rihaio ~ 
160( 1 1 ^) 3 cm- 3 (A/0.05)- 3 in a disk (where A is the spin parameter). A supernova remnant at high 
z expands into a denser ISM: since it spends a shorter time in the Taylor-Sedov phase, most of 
its energy is radiated at a smaller radius, where the effective temperature is higher. This implies 
a harder spectrum for thermal emission. In addition, the density and temperature of gas in star 
forming regions is determined by the properties of H2 cooling (which saturates at n ~ 10 4 cm" 3 and 
T ~ 300K), rather than metal cooling as in the local universe, (hi) Potential wells are significantly 
shallower, so pressurised regions (hot gas, strong magnetic fields) cannot be efficiently confined, (iv) 
The CMB energy density Ucmb oc (1 + z) 4 , so inverse Compton radiation becomes particularly 
efficient at high redshift. Because of (hi) and (iv), relativistic electrons cool predominantly by 
inverse Compton scattering rather than synchrotron emission. Since inverse Compton emission is 
likely to be the most promising mechanism for X-ray emission, I shall consider it at length. 



2.3. Inverse Compton emission 

The energy loss rate for a relativistic electron with Lorentz factor 7 is given by: 

Vic = \o TC1 2 U rad = 1.12 x 10- 16 (^) 4 (t^) 2 ergs- 1 (4) 

for U ra d = JTcmb- The loss rate by synchrotron radiation is given by substiting Ub for U ra d, 
^synch = \^tc^ 2 Ub- In the local universe, galaxies with relatively quiescent star formation emit 
most of their electron energy in synchrotron radiation. Starbursts in the local universe can radiate 
efficiently in IC, as the energy density in the local radiation field is sufficiently high (typically, 
U r ~ 10 -8 erg cm -3 as opposed to U r ~ 10~ 12 ergcm" 3 in our Galaxy). Seed photons are provided 
by IR emission from dust grains. By contrast, at high redshift, all star forming regions will emit 
in inverse Compton radiation, as the CMB provides a universal soft photon bath of high energy 
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density, Ucmb = 4 x lCT 9 (ij^) 4 ergcm~ 3 . An electron with Lorentz factor 7 will boost a CMB 
photon of frequency v to a frequency v = 7 2 f Q . Thus, the rest frame frequency of a CMB photon 
(at the peak of the blackbody spectrum) which undergoes inverse Compton scattering is: 

Note that the observed frequency is independent of source redshift, since the higher initial frequency 
and redshifting effects cancel out. 

Below, I examine in detail the mechanisms by which relativistic electrons lose energy, to see if 
indeed inverse Compton radiation will predominate. 



2.3.1. Energy loss processes for relativistic electrons 

Once relativistic electrons are produced, they can cool via a variety of mechanisms. Let us 
examine them in turn (for more details see Pacholczyk 1970, Daly 1992). 

Synchrotron radiation The relative emission rate in inverse Compton and synchrotron 
emission is given simply by the relative energy densities in the radiation and magnetic fields, 
EicjEsyn = U rad /U m . This yields: 

j££ = L1 x (syVinV (6) 

Note that in energy loss terms the CMB may be characterized as having an effective magnetic field 
strength BcMB,eff = 3.24 x 10 2 (^f ) [J.G. Could magnetic fields in proto-galaxies possibly reach 
these high values? A reasonal assumption is that Pb ~ Prei < Pgas, where Pb is the magnetic 
field pressure, P re i is the pressure in relativistic particles, and P gas is the thermal gas pressure. 
Local observations of synchrotron and inverse Compton emission from radio galaxies are consistent 
with equipartition Pb ~ P re i (Kaneda et al 1995). I have assumed that the energy injection into 
relativistic particles is ~ 10% of the total kinetic energy of a supernova, so P re \ < P gas should be 
a strict upper bound. Thus, Pb < Pgas gives the upper bound: 

B < 6 (l^) 1/2 (lo^)^ < 7 > 

where n is the baryon number density (note that gas with temperatures > 10 4_5 K will escape 
from the shallow potential wells of the first proto-galaxies). If the magnetic field exceeds the above 
value, the over-pressurised lobe will expand on the dynamical time scale until the magnetic pressure 
drops. Thus, for z > 5, it seems likely that synchrotron energy losses will be unimportant. 

Ionization & Cherenkov losses Interactions with the non-relativistic gas will result in 
energy losses via ionization and Cherenkov emission of plasma waves at a rate independent of the 
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Lorentz factor, Ei on ps 9 x 10 19 nergs 1 , which implies that the relative energy loss rate is: 

/ n \-i/1-i-^\ 4 /^'\2 

lic^Fs; v 



(i^)" I (^) 4 (^) : 



n \i/2 / l + z w2 



Thus, for 

7 > 76reafc ~ 100 ( — ^) ' (^— -)" 2 (9) 

Vlcm ci / 10 

inverse Compton losses are more important than ionization losses. The cutoff Lorenz factor at the 
lower end 7 co , is of interest since the lower end accounts for most of the electrons, both in terms of 
number and energy: N e i ectrons (> 7) oc 7 _2a = 7" L6 ; E electrons (> 7) oc <y- 2a+1 = 7~ - 6 (q = 0.8). 
Note that the electrons that cool via ionization losses do not drop out completely, but merely form 
a flattened distribution with 7^ ~ 73 — 350n(t/10 yr). These electrons with lower Lorentz factors 
could scatter CMB photons to optical and UV frequencies. 

Free-free radiation Free-free radiation results in an energy loss rate Ef ree _f ree = 6 x 
10 -22 ri7ergs _1 . Thus, free- free radiation only dominates over ionization and Cherenkov radation 
for 7 > 1500. However, in this regime inverse Compton losses dominate, since 

^ = 190 M^r 1 f i±^y m (io) 

Thus, free-free emission is never important in cooling relativistic electrons at high redshift. 

In summary, the inverse Compton emission is the dominant energy loss mechanism between a 
lower and upper energy cutoff. From equation (9), the lower frequency break is determined by the 
competition between the inverse Compton loss rate and ionization and atomic cooling losses, below 

Elower = JbreakhvCMB = 70 ( ^-3 ) (~Yo~) gV ' ( U ) 

where I assume the seed photon ucmb = 1.6 x 10 12 (^p ) GHz lies at the peak of the CMB blackbody 
spectrum. The competition between inverse Compton cooling losses and the rate of energy injection 
by Fermi acceleration determines the upper energy cutoff. The timescale for losses by inverse 
Compton radiation is: 

tlj/e = § = 7.9 x 10* (I^^i)" 1 years (12) 

Equating the Fermi acceleration timescale t acc ~ r L c/v 2 sh = 1.3x 10~ 3 (ggg) (j^g) Qoookms-O ^ 
(where is the Larmour gyroradius, and v s h is the typical shock velocity) to tnf e (equation (12), 

I obtain for the maximum Lorentz factor j ma x = 7.4 x 10 6 ^j^q^ ^ ( 2000km s- 1 ) (tct) 2 wmcn 
corresponds to an upper energy cutoff: 



E „„_ 360 (^)(_ S _) 2 (i±£)- 3 GeV . (l 3) 
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In section (2.3.2), I derive the form of the spectrum. For the power law spectrum obtained, 
L u oc v _1 , the specific luminosity depends only logarithmically on the energy cutoffs: L u = 

LtQt y~ 1 

log(^u P por/^lowcr) 



2.3.2. Inverse Compton spectrum 

The Fermi shock acceleration mechanism for cosmic rays involves a steady growth of particle 
energy as a particle scatters back and forth across the shock front. This naturally produces a power 
law electron energy spectrum with energy spectrum dn/d'j oc 7" p , where p = (x + 2)/(% — 1) and 
X is the compression ratio for the shock (e.g., Jones & Ellison 1991). The final emission spectrum 
L v oc v~ a depends on the electron energy spectrum; the synchrotron or inverse Compton spectral 
index is a = (p — l)/2. Thus, the spectral index is determined by the shock structure rather than 
the details of the scattering process. For a strong adiabatic shock \ = 4 and dn/d'j oc 7~ 2 ; for an 
isothermal shock, \ ^ 1 an d dn/d'j oc 7 -1 . The assumption of adiabaticity is most appropriate for 
the Taylor-Sedov phase, when most of the electrons are accelerated. Galactic SNR show a mean 
radio spectral index a = 0.5 ± 0.15 (Droge et al 1987), which agrees with dn/d'j oc 7 -2 . I shall 
use this as the canonical electron injection spectrum in this paper. The steepening of the diffuse 
synchrotron emission to a ~ 0.8 is likely to be due to energy losses as the electrons age. I examine 
the emission spectrum in detail below. I assume that there is one supernova for every 100 M & of 
stars formed, that each supernova liberates ~ 10 51 erg in kinetic energy, and ~ IMq of metals. 

The equation for the evolution of the electron population is given by (Ginzburg & Syrovatskii 
1964, Sarazin 1999): 

= ^[6(7)^(7)1+0(7) (14) 

where Ni^j) is the number of electrons in the range 7 to 7 + dj, the rate of production of new 
relavistic electrons is given by Q{j), and the rate of energy loss of an individual particle is given 
by 6(7) = d'j/dt. From equation (12), we see that the electron lifetime is shorter than both the 
Hubble time for the redshifts under consideration z < 30 and typical timescales for starbursts 
(~ 10 7 years). We can thus assume that the electron population quickly reaches a steady state 
where the energy losses due to inverse Compton scattering (for high 7) and Coulomb collisions (for 
low 7) balance the injection rate due to supernova explosions, and set the time derivative to zero. 
Since injected electrons survive for less than a Hubble time, I ignore evolution in loss rate due to 
the evolution of the CMB energy density. Assuming that each supernova injects a population of 
electrons with ^(7) = Qo7~ p , where p ~ 2, then the steady state solution to equation (14) is given 
by: 

Nb) = 4xl „ 61 (l±i)- 4 (_gL I ) T -( P+1 , ;7>1W (15) 
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where ^break (equation (9)) is the transition energy between the regimes where ionization and inverse 
Compton losses dominate. Thus, the electron population flattens by one power at low energies and 
steepens by one power at high energies, as compared with the distribution function for the injected 
population. Since a = — (p — l)/2, the emitted spectrum flattens by a half power and steepens by 
a half power at low and high energies respectively. In particular, for p = 2, as is appropriate for 
adiabatic shocks, L u oc const at low energies and L u oc u^ 1 at high energies. 

The spectrum of inverse Compton radiation is given by (Sarazin 1999): 

L v = 127ra T f 7V( 7 )d 7 / J{-^ r )F{x)dx (16) 
Ji Jo 4rfx 

where N(j) is the electron energy spectrum, J[y) = B u (Tcmb,o(^ + z )) is the CMB blackbody 
spectrum, and F(x) = 1 + x + 2xln(x) — 2x 2 . This differs only by a normalization correction from 
assuming that all photons reside at the peak of the blackbody spectrum, v = 1.6 x 10 11 (1 + z)Hz. 
This is because the seed photon spectrum is narrow compared to the electron energy distribution. In 
figure (1), I show example spectra computed with equations (15) and (16). Note that the efficiency 
of Coulomb cooling depends on the assumed electron density; I assume that n e ~ Sn (l + z) 3 , 
where I assume the overdensity to be either 5 ~ 200 (to mimic the overdensity at the virial radius) 
or 5 ~ 10 4 (to mimic the overdensity for collapsed gas in dense star forming regions). At z ~ 10, 
this assumption corresponds to n e ~ 4 x 10~ 2 cm -3 and n e ~ 2 cm -3 respectively. Also plotted for 
reference is the spectrum of the same starburst with a Salpeter IMF, assuming an escape fraction 
of ~ 1%, and the spectrum of a mini-quasar with spectrum L u oc z^~ 18 , normalized to have the 
same energy release above the Lyman limit as the starburst. Note that the inverse Compton case 
has the hardest spectrum of all. Since the spectrum is so hard, photoelectric absorption by the 
host galaxy does not significantly attenuate the ionizing flux. Most UV photons produced will 
not escape the host galaxy, whereas X-ray photons with E > 270 ( w -2i"^--2 ) eV where Nhi is the 
column density in the host galaxy, will escape unimpeded. Since uL u ~ const, most of the energy 
in ionizing photons escape. Note, incidentally, that the escape fraction for UV photons produced 
by inverse Compton may be significantly higher than photons of the same frequencies produced 
by stars, as relativistic electrons can disperse from the star forming region, where gas densities are 
highest and most photon captures take place. 



2.4. Zero-metallicity star formation 

The stellar IMF at high redshift under conditions of low or zero metallicity is unknown. Up 
to now, I have assumed a Salpeter IMF, in which one supernova expodes for every ~ 100 M Q of 
stars formed, and each supernova deposits on average ~ 10 50 (g^-)erg in relativistic electrons and 
~ 1 Mq of metals. Since the estimated amount of star formation at high redshift is calibrated to 
the observed IGM metallicity at z ~ 3, it is worth asking whether the energy injected into X-rays 
per solar mass of metals produced, ez ~ 10 50 ergMQ 1 , could change significantly at high redshift. 
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I have argued that variations in the IMF will not effect significant changes in ez or the ratio of 
energy emitted in stellar UV ionizing photons to inverse-Compton X-rays, since the same massive 
OB stars which produce UV ionizing photons also explode as supernova, producing both relativistic 
electrons and metals. 

However, zero-metallicity star formation may result in very different stellar populations from 
that seen at low redshift. The lack of efficient cooling mechanisms could result in extremely top 
heavy stellar IMFs (Larson 1998, Larson 1999) and in particular the production of "Very Massive 
Objects" (VMOs) in the range 10 2 - 10 5 M (Carr, Bond & Arnett 1984). Such stars have high 
effective temperatures and produce a much harder spectrum than ordinary stars (Tumlinson Sz 
Shull, 2000, Bromm et al, 2000). Moreover, the distribution of stellar endpoints is very different 
from a normal IMF (Heger, Woosley & Waters 2000). Stars with masses between 10 and 35 
Mq explode as type II supernovae. While it is not known whether metal free stars with masses 
between 35 and 100 M Q will explode-they might collapse to form black holes-stars with masses 
~ 100 — 25OM are disrupted by the pair production instability, once again producing an energetic 
supernova event and dispersing metals. Stars more massive than 25OM should collapse completely 
to black holes, without ejecting any metals (unless they eject their envelopes during hydrogen shell 
burning, in which case it is possible for them to explode). The latter provides an obvious mechanism 
for seeding supermassive black holes to form AGN. A number of studies of zero-metallicity star 
formation suggest that the Jeans/Bonner-Ebert mass and hence the lower mass cutoff for star 
formation is very high, M* > 10 2 - 10 3 M & (Abel, Bryan & Norman 1999, Padoan, Nordlund & 
Jones 1997). 

There are three main observations to make: (i) VMOs with initial stellar masses 100M Q < 
M < 25OM which are disrupted by the pair-production stability show ez(PP) ~ e^(type II), 
which imply that even if these objects are abundant in the early universe, our estimates of the 
level of inverse Compton X-ray emission do not change strongly. In particular, -Explosion = 
6.3 x 10 52 ( io^m ) mm[0-i3,(i^ 4>l) — 1 er g w ]^^ e ^he yield in elements heavier than helium is Z c j = 
min [(1 — <Pl), 0.5], where (pi is the fraction of the initial mass lost during hydrogen burning 
(Carr, Bond & Arnett, 1984). Thus, e z (PP) = 1.3 x 10 50 [j^j^) * ' § HlM^M^Mif 1 erg Mq 1 « 
e z(type II). (ii) Stars which directly collapse to form black holes represent an additional, unac- 
counted source of ionizing photons, since they are not included in the metal pollution budget (of 
course, their most important contribution could lie in seeding AGN formation), (iii) Zero metal- 
licity is a singularity: true zero-metallicity stellar populations differ greatly from low metallicity 
ones. Even the introduction of trace amounts of metals Z ~ W~ 4 Zq introduces important changes 
in stellar structure and evolution (Heger, Woosley & Waters 2000). Furthermore, trace amounts 
of metals drastically reduces the abundance of VMOs by allowing efficient cooling past the 300K 
barrier imposed by H2 cooling, reducing the Jeans mass and thus the minimum mass for star for- 
mation. Prompt initial enrichment is quite plausible: the lifetime of massive stars before they 
explode to pollute the IGM with metals is ~ 10 6 years, which is a short fraction of the Hubble time 
t//~8x 10 8 (i^)~ 3 / 2 yrs even at high redshift. Although the degree of metal mixing is uncertain 
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(e.g. Gnedin & Ostriker 1997, Nath & Trentham 1997, Ferrara, Pettini & Shchekinov 2000), note 
that first star forming regions are very highly biased, and subsequent generations of the halo hi- 
erarchy collapse in proto-clusters and filaments very close to the first star clusters. Thus, stars of 
finite metallicity could quickly predominate, even if the mean metallicity of the universe is close to 
primordial (Cen & Ostriker 1999). It is therefore possible and even likely that true zero metallicity 
star formation was confined to a very small fraction of the stars formed at high redshift, and thus 
negligible in terms of the energy budget for reionization. 



3. Observational signatures 

In this section, to estimate number counts I use the Press-Schechter based high redshift star 
formation models of Haiman & Loeb (1997), which are normalised to the observed metallicity of 
the z=3 IGM, Z = 10~ 3 — 1O~ 2 Z . In these models, in every halo capable of atomic cooling (i.e., 
with a virial temperature T v i r > 10 4 K), a fixed fraction of the gas f s t ar = 1.7,17% (for Z(z = 
3) = 10~ 3 , 1O~ 2 Z respectively) fragments in a starburst lasting t Q ~ 10 7 yrs. The correspondence 
between halo mass and star formation rate is thus SFR rs 2 ( io«Mq ) (tTIt) M©yr _1 , and each 
halo is only visible for some fraction of the time *? . 



3.1. Hell recombination lines 

One possible signature of the hard spectrum produced by inverse Compton X-rays would be 
Hell recombination lines from the host galaxy. Indeed, such lines may well be detectable from 
the first luminous objects with sufficiently hard spectra, such as mini-quasars or metal- free stars 
(Oh, Haiman & Rees 2000). The different sources may perhaps be distinguished on the basis 
of line widths and line ratios (Tumlinson, Giroux & Shull, 2000). However, inverse Compton 
emission produces too few Hell ionizing photons for such recombination lines to be detectable with 
NGST at high redshift. For a Salpeter IMF, N ion ,Hi = ^( im™- 1 ) photons s" 1 from stellar 
UV radiation. On the other hand, since secondary ionizations of Hell are negligible (Shull & 
van Steenberg 1985), the production rate of Hell ionizing photons from non-thermal emission 
is Nion,Heii = f££ n j£ ~ 2.5 x 10 49 ( 1A fJ^-i ) photons" 1 , where v thin is the frequency at 

which the halo becomes optically thin. Thus, Q = N^^/N^ ~ 2.5 x 10~ 4 , as compared with 
Q ~ 0.05 for stellar emission from metal-free stellar population, and Q ~ 4r a = 0.08 — 0.25 for 
a QSO, where L v oc v~ a and a = 1 — 1.8, implying that the relative flux in Hell and Ha, (5 
recombination lines is much smaller for inverse Compton emission than metal-free stars or AGN. 
The luminosity in a helium recombination line i may be estimated as L\ = Qf\L}i a , where _Lh« 
is the Balmer a line luminosity, f\ = -2*— [ a ?}i i J}\ lHI1 — ), e is the fraction of SN energy which 

J ' J1 ]Hu \a B (HeII)n He iii ) ' toJ 

emerges as IC emission, and ji,jna may be obtained from Seaton (1978). The observed flux is 
Ji = ^tkl = O-O 4 (<&) (TOT' ((FT) (TOT)) (li^) nJ y where R is the s P ectral resolution, 
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and qi = fcvua/vi, where / A 4686 = 0.74, / A1640 = 4.7, A3203 = 0.30. This is undetectable by 
NGST, which in 10 5 s integration time requires for a 10 a detection a flux .F(10<7) ~ 30nJy at these 
frequencies (observed wavelengths 1 < A < 5.5//m) and spectral resolution R ~ 1000 (see Oh, 
Haiman & Rees 2000 for details). 

It is worth mentioning that a hard source may produce relatively little recombination line 
flux and yet play an important role in reionizing the universe. The recombination line flux from 
a source is oc (1 — f esc ), whereas the ionizing flux escaping into the IGM is oc f esc (where f esc 
is significantly higher for hard sources: all photons with E > Eh a io,thin = 270 ( 10^^-2 ) e ^ can 
escape freely from the host). Secondary ionizations are unimportant within a host halo: much of 
the gas in fully ionized, in which case the energetic photoelectron deposits its energy as heat. In any 
case, the halo is optically thin to the most energetic photons E > Eh a i jhin which are important for 
secondary ionizations. Thus, in recombination line flux what matters is the total number of ionizing 
photons produced, whereas for the reionization of the IGM what matters is the total output energy 
(since in a largely neutral medium with x e < 0.1 the total number of ionizations for a photon of 
energy E p h oton , including secondary ionizations, is iVj on ~ Eph oton /37eV (Shull & van Steenberg 
1985). Note, however, that as the medium becomes more ionised an increasing fraction of the 
energy is deposited as heat, and an additional source of soft photons is needed for reionization 
to proceed (Oh 2000a)). Thus, for instance, for f esc ~ 1%, inverse Compton radiation makes a 
negligible contribution to the Ha luminosity, s "u a r ~ < T st°iiar = 2-5 x 10 , but the ionizing IC 
and stellar radiation escaping from the host are energetically comparable (from equations (1) and 
(3), L^p « Lfjy lar ), and thus they produce roughly equal number of ionizations in the IGM. 



3.2. X-rays and gamma-rays 

Unfortunately, direct detection of inverse Compton X-rays from high-redshift star clusters is 
unlikely. From equation (1), the observed flux in X-rays from a star forming region is: 

f = 7T = 5 x 10 2U — — T erg s 1 cm 1 (17) 

J Audi vo.i/ V 1M ©y r / V 10 y v ; 

By contrast, the Chandra X-ray Observatory (CXO) sensitivity (see http://chandra.harvard.edu/) 

for a 5a detection of a point source in an integration time of 10 5 s in the 0.2-10 keV range is Fx = 4x 

10~ 15 erg s _1 cm -2 . The next generation X-ray telescope, Constellation-X (http: //constellation. gsfc.nasa.gov/) 

is optimised for X-ray spectroscopy and will not have significantly greater point source sensitivity: 

for a similar bandpass and integration time it will be able to detect objects out to a flux limit of 

Fx = 2 x 10~ 15 ergs -1 cm~ 2 . Thus, in a week CXO and Constellation-X will at best be able to 

detect starbursts (with seed photons provided by the IR radiation field from dusty star forming 

regions) with star formation rates SFR ~ 100M Q yr _1 (corresponding to Lx ~ 4.8 x 10 42 ergs _1 ) 

out to redshift z ~ 1 (note that the Lyman-break galaxies detected at z ~ 3 (Steidel et al 1996) 

typically have have inferred star formation rates SFR ~ 100M Q yr _1 ), while the very brightest 
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starbursts, with star formation rates of SFR~ lOOOM yr~ , might be detectable out to z ~ 3. As 
discussed in section 3.4, simultaneous detection of synchrotron radiation with the Square Kilometer 
Array will then constrain magnetic field strengths in these objects. 

Detection in gamma ray emission is also unlikely. The upcoming Gamma-ray Large Area 
Space Telescopy (GLAST) (http:/ /glastproject. gsfc.nasa.gov/) will have a flux sensitivity of ~ 2 x 
10~ 9 photons cm~ 2 s _1 in the 20MeV-300GeV range, whereas star-forming regions at high redshift 
will have fluxes of at most K, ~ 10~ 15 (^y) ( 1M ^ F -^-i ) (^) 2 photons cm^s -1 . 

What fraction of the unresolved X-ray and gamma-ray background could be due to star for- 
mation at high redshift z > 3? Observations of absorption in CIV and other metals in Lya forest 
absorption lines with 10 14 - 7 cm _2 < N HI < 10 16 cuT 2 indicate that Z ~ lCr 2 Z ~ 2 x 10~ 4 . Each 
supernova produces about ~ 1M Q of metals (Woosley & Weaver 1995). Assuming (ftb,h) = 
(0.04,0.65), this implies one supernova every ~ 1000 kpc 3 (comoving). If each supernova in- 
jects ~ 10 50 ergs in hard X-rays, the comoving X-ray energy density was Ux ~ (l+s) ~ 
2 x 10 _6 eVcm -3 . If the mean source redshift was z ~ 5, the distant sources produced a dif- 
fuse flux of J ~ lkeVs _1 cm _2 sr _1 . One can make a more detailed estimate using the equation of 
cosmological transfer (Peebles 1993): 

where z Q is the observer redshift, v = u (l+z)/(l + z ), and e(v, z) is the comoving X-ray emissivity. 
Using the star formation model of Haiman & Loeb (1997), I obtain the estimate: 

vJ v = 0.8 keVs- 1 cm- 2 sr- 1 for z source > 3 (19) 

In Fig (2), I display this predicted level of emission against the observed X-ray and gamma-ray 
background in the 1 keV -100 GeV range, from analytic fits to the ASCA and HEAO A2,A4 data 
in the 3-60 keV range (Boldt 1987), the HEAO 1 A-4 data in the 80-400 keV range (Kinzer et al 
1997), the COMPTEL data in the 800 keV - 30MeV range (Kappadath et al 1996). Data points 
from EGRET in the 30 MeV - 100 GeV range (Sreekumar et al 1998) are also shown. The level of 
the unresolved background of course depends on the sensitivity of the instrument; recently CXO 
resolved ~ 80% of the hard X-ray background in the 2-10 keV range into point sources (Mushotzky 
et al 2000). This agrees well with the predictions of XRB synthesis models (e.g., Madau, Ghisellini 
& Fabian 1994), which use AGN unification schemes to reproduce the observed spectral shape 
of the XRB. A prediction of the IC scenario presented here is that a non-trivial fraction of the 
X-ray /gamma-ray background will not be resolved into point sources with upcoming missions, due 
to the extreme faintness of high redshift sources. 

It is particularly intriguing that both the amplitude and spectral shape of the gamma-ray 
background as observed by EGRET is well-matched by the predicted level of gamma-ray emission 
in this model. This raises the exciting possibility that the majority of the observed gamma-ray 
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background comes from inverse Compton emission at high redshifts. At present, the origin of the 
gamma-ray background is still unknown. The most favoured scenario for some time was that it 
is due to unresolved gamma-ray blazars (Bignami et al 1979, Kazanas & Protheroe 1983, Stecker 
& Salamon 1996): the observed blazar 7-ray spectrum has an average spectral index compatible 
with the observed GRB (Chiang & Mukherjee 1998). However, extrapolation of the observed 
EGERT blazar luminosity function implies that unresolved blazars can account for at most ~ 25% 
of the diffuse 7-ray background (Chiang & Mukherjee 1998). The unresolved blazar model will 
most likely be decisively tested by GLAST (Stecker & Salamon 1999), which will be two orders 
of magnitude more sensitive than EGRET. A host of other models include pulsars expelled into 
the halo by asymmetric supernova explosions (Dixon et al 1998, Hartmann 1995), primordial black 
hole evaporation (Page & Hawking 1976), supermassive black holes at very high redshift (Gnedin 
& Ostriker 1992), annihilation of weakly interactive big bang remnants (Silk & Srednicki 1984, 
Rudaz & Stecker 1991) and finally, inverse Compton radiation from cosmic ray electrons in our 
own Galaxy (Strong & Moskalenko 1998, Dar & De Rujula 2000), and from collapsing clusters 
(Loeb & Waxman 2000). However, to date the possibility of inverse Compton emission from high 
redshift super novae has not been discussed. 

A scenario in which the majority of the gamma-ray background comes from inverse Compton 
emission at high redshift make a number of firm predictions: (i) As previously mentioned, the 
majority of the GRB will remain unresolved by GLAST, due to the extreme faintness of the 
contributing sources, (ii) After removal of the Galactic contribution, which is correlated with 
the structure of our Galaxy and our position within it (Dixon et al 1998, Dar et al 1999), a 
highly isotropic component of the GRB will still be present, (iii) The GRB should be extremely 
smooth, and exhibit significant fluctuations only at extremely small angular scales. The fluctuations 
should be dominated by the Poisson rather than the clustering contribution (see Oh 1999). Using 
(S n ) = f Q 6c ^S n , where S c is the cut-off flux for point source removal, and the star formation 
model of Haiman & Loeb (1997), I find that for z > 3 sources (which are too faint to be removed as 
point sources), — = 3.4 x 10~ 2 (J^) \ where I is surface brightness (the angular resolution of 
GLAST is expected to be of order 1 — 5 arcmin). (iv) The gamma-ray background at E > lOOGeV 
should be attenuated, due to pair production opacity against IR/UV photons (e.g., Salamon & 
Stecker 1998, Oh 2000b). High energy photons initiate an electromagnetic cascade which transfers 
energy from high energy photons to the lower energy portion of the spectrum, where the universe is 
optically thin (Coppi & Aharonian 1997). Note that the EGRET spectrum was directly determined 
with data only up to 10 GeV; beyond 10 GeV larger uncertainties exist due to backsplash in the 
Nal calorimeter, and Monte-Carlo simulations were used to determine the differential flux in the 
10-30,30-50 and 50-120 GeV range (Sreekumar et al 1998). Thus, the EGRET data points with 
E > lOGeV are less reliable. It would be intriguing to see if GLAST (sensitive out to 300 GeV) 
indeed shows an absorption edge to the gamma-ray background at higher energies. It would also be 
interesting to look for absorption of the gamma-ray background in a line of sight passing through 
a massive cluster (indicating that the gamma-rays come from higher redshifts than the cluster); 
the pair-production optical depth through a cluster is of order r ~ 2n 1 aTr v i r ~ 0.4 (assuming 
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0.1cm 3 and r v i r ~ IMpc). 



3.3. 



CMB constraints 



Will the upscattering of CMB photons by relativistic electrons at high redshift produce an 
observable signal in the CMB, or violate any present observational constraints? When the electrons 
are relativistic, CMB photons are inverse Compton scattered to such high energies (completely out 
of the detector bandpass) the process may simply be thought of as absorption, 51 /I ~ rj e/ , where 
T r e el is the optical depth of relativistic electrons, and / is the number flux of CMB photons. The 
flux decrement due to the absorption of CMB photons may be estimated as: 



where Jcmb is the blackbody surface brightness of the CMB, and iVJ is the total number of 
relativistic electrons in the system. The steady state number of relativistic electrons is given by 
N£ el = f djN(j) where N(j) is given by equation (15). Since the CMB flux peaks at v ~ 10 n Hz, 
the absolute magnitude of the flux decrement is maximized by going to similar frequencies. At 20 
GHz, the highest frequency detectable by SKA, the SKA has an rms sensitivity of ~ 6nJy for a 10 5 s 
integration. On the other hand, the flux decrement is AS„(20GHz) ~ 3.5 x 10~ 5 (±±f-) 2 (N™ 1 /3 x 
10 58 electrons) nJy, which is unobservably small. 

One might hope to detect the mean signal of all the CMB photons upscattered by relativistic 
electrons. In particular, since the number of CMB photons is no longer conserved, the absorption 
might be detectable as a chemical potential distortion of the CMB. Let us estimate the number of 
CMB photons destroyed. Each supernova upscatters at most N scatte red ~ tE$N / (Ex) ~ 10 60 CMB 
photons, where I have set the average photon energy (Ex) ~ lOOeV (note that since the number 
of photons N u oc v^ 1 , most of the upscattered photons are of low energy). For a metallicity of 
Z ~ 10~ 2 Z & at z ~ 3, one supernova has gone off every comoving Vsn ~ 1000 kpc 3 , and thus the 
comoving number density of upscattered CMB photons is 5n ~ N scattere ^/VsN ~ 4 x 10 _8 cm 3 . 
Since n 7 ~ 400cm -3 , we have 5n/n ~ 10~ 10 , which results in an undetectably small chemical 
potential distortion. Thus, the upscattering of CMB photons at high redshift does not violate any 
distortion constraints on the CMB. 

If the IGM is reionized inhomogeneously, as in canonical models, then secondary CMB anistropies 
will be created by CMB photons Thompson scattering off moving ionized patches (Agahanim et al 
1996, Grusinov & Hu 1998, Knox et al 1998). The power spectrum is generally white noise, with 
AT/T ~ 10 -6 — 10~ 7 , peaking at arc-minute to sub arc-minute scales. However, if the IGM is reion- 
ized fairly homogeneously by X-rays then over a line of sight the positive and negative contributions 
of the velocity field will cancel out. In this case, only the second-order Ostriker-Vishniac effect due 
to coupling between density and velocity fields will be present. A null detection of the inhomoge- 
neous reionization anisotropy could place upper limits on the patchiness of reionization, although 




'A 
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this will be a difficult measurement as the inhomogeneous reionization and Ostriker-Vishniac signals 
are likely to be of comparable strength (Haiman & Knox 1999). 



3.4. Radio observations 

The supernovae that exploded generate magnetic turbulence and magnetic fields, which allow 
Fermi acceleration to take place (Jones & Ellison 1991, Blandford & Eichler 1987); observations of 
local radio galaxies in non-thermal radio and X-ray emission yield field strengths consistent with 
equipartition between relativistic particles and the magnetic field (Kaneda et al 1995). If such 
magnetic fields are present in the first star clusters, the relativistic electron population is a source 
of synchrotron radio emission as well as inverse-Compton emission. Below I find that for a given 
relativistic electron population, the Square Kilometer Array (SKA) will be much more sensitive 
to non-thermal radio emission than CXO or Constellation-X will be to inverse Compton X-ray 
emission. Radio observations will thus allow one to establish the presence of relativistic electrons 
in objects too faint to observe directly in X-ray emission. Since one knows the CMB energy density 
exactly as a function of redshift, given reasonable assumptions for the magnetic field the observed 
radio emission allows one to immediately estimate the amount of inverse Compton X-ray emission 
which must be taking place, 

t _ Ucmb(z) t , . 

l^X — 77 ^synch (/1J 

The transition Lorentz factor jbreak & t which electron energy losses are dominated by inverse 
Compton rather than ionization losses is given by equation (9). Above an observed radio frequency 

"break = "Ll 2 break /(1 + *) = 2 - 8 * lO 4 ^) (uZF*) (tuT* Hz ( where v l = ^ is the electron 
gyrofrequency) , the steady state electron population is determined by balance between supernova 
injection and inverse Compton losses, and is given by equation (15). From standard formulae for 
synchrotron emissivity (Rybicki & Lightman 1979) ef^ nc ^ = crj^ 2 (3 2 U mag cn(E)dE/di; (where n(E) 
is the number density of emitting electrons in dE) this yields a synchrotron luminosity: 

LT < _ , 7 „ lfl » (_b_) - ( _y - ( £ } /(l> Vi) ergs - lHz - 1 (22) 



where f(z, f/ 7 ) = min 



(the latter term in brackets is used if the 



V 10 J ' \^4.2xl0- 9 crg cm- 3 J 

stellar radiation field has a higher energy density than the CMB). By contrast, the thermal free-free 
emission is given by (Oh 1999): 



(23) 



Thus, assuming a = 1, at observed frequencies v < u trans = 2.3 (^jjp) 1 ^i^g) (or) f( z ' ^r) 
synchrotron emission dominates over free- free emission. This is well within the 0.1-20 GHz capabil- 
ity of the SKA. However, this is only true if the power law for electron population extends to high 
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energies. If it is truncated at some maximum Lorentz factor ^ max , then synchrotron emission is 

only observable for v < u max = VLl max = °- 28 (tt) 1 (t^g) i 1 ^) 2 GHz - If v rnax < vtrans, 
this will be manifested by an abrupt drop in radio flux at v max , beyond which the emission 
takes the flat spectrum free-free emission form (this does not take place for expected values of 

1 /2 

Imax = 7.4 x IO^y^-) ( 2000^18-0 (tct) " ^ see ec l uatio n ( 13 ))- Observation of such a drop 
will yield valuable constraints on B, "y ma x in the first objects; the value of ^ ma x in turn constrains the 
upper energy cutoff in the inverse Compton X-ray /gamma-ray spectrum, which could be important 
in determining whether high-redshift starbursts make a significant contribution to the gamma-ray 
background. 

Can radio emission from star clusters at high redshift be detected by the proposed Square 
Kilometer Array? Non-thermal emission can be distinguished from free-free emission with multi- 
frequency observations to identify frequency regimes where the spectral slope is steep. The SKA 
detector noise may be estimated as: 

s — = T^ms - 25 (ass*) (ik) < 24 > 

where I have used A c g/T sys = 2 x 10 8 cm 2 /K for the SKA (Braun et al 1998), and I assume a 
bandwidth Au « 0.5za The flux density due to non-thermal emission from a high-redshift star 
cluster, assuming a = 1, is: 

s, - ^^--^r (i^) (£) m~ 2 n^ 

(25) 

Thus, a source with SFR~ lOMoyr" 1 at z ~ 9 can be detected as a 10<r detection in 10 days. 
To estimate the number of sources detectable by SKA, I use the Press-Schechter based high red- 
shift star formation models of Haiman & Loeb (1997), and define an efficiency factor f ra dio = 

( lO^G ) (ot) ({frf)' where f s t ar is the fraction of halo gas which fragments to form stars. In 
Figure (3), I display the number of sources above a given redshift which may be detected in non- 
thermal emission in the 1° SKA field of view, assuming f ra dio = lj 0.01. Also shown is the number 
of sources which can be detected in free-free emission at 4 GHz. One should be able to detect a 
large number of sources at high redshift, z > 5. Thus, radio observations of non-thermal emission 
can serve as a useful proxy for X-ray observations in allowing one to estimate Lx, and thus the 
overal level of X-ray emission at high redshift. 



3.5. Multi-wavelength observations 

While a detailed study of high-redshift multi-wavelength campaigns is beyond the scope of this 
paper, below I describe some possible follow-up observations if synchrotron emission is detected at 
high redshift. 



-22- 



• Redshift estimates Thus far the most efficient way to select high-redshift radio galaxies 
has proven to be the observation of steep radio spectra a < —1.3 (Chambers et al 1990, van 
Breugel et al 1999). This is at least partially due to the fact that extremely bright radio 
galaxies have SEDs which steepen with frequency; the k-correction then implies that sources 
at increasing redshift have steeper observed spectra. Steepening with redshift due to inverse 
Compton losses, as well as selection effects (i.e., brighter sources have stronger magnetic fields, 
and thus more rapid synchrotron losses) could also play a role (Krolik & Chen 1991). This 
technique may fail for fainter sources at high redshift since for v > min(fj raras , v m ax\ the radio 
flux will be dominated by free-free emission and the spectra will appear flat. For instance, 
for B < 3 (^j(f) 2 5 (ot)° 5 A*G, we have vtrans < 160MHz and SKA will only detect free-free 
emission within its frequency coverage. An efficient way to select high-redshift objects prior 
to reionization would be to perform broad-band deep field imaging with NGST and select 
Lyman-break dropouts as has been done at z ~ 3 (Steidel et al 1996); before the epoch of 
reionization one must select 'Gunn- Peterson dropouts', i.e. galaxies with no flux shortward of 
rest-frame HI Lya. Yet another method of selecting high-redshift objects would be to perform 
a joint survey in the submm with the Atacama Large Millimeter Array (ALMA); at low 
redshift the submm dust emission scales almost linearly with other star formation indicators, 
such as radio and UV emission. However, for z ~ 0.5 — 10 and a given dust emission SED the 
K-correction almost balances the cosmological dimming of a source, implying a flux density 
almost independent of redshift (e.g., Blain et al 2000). This has spawned suggestions to 
obtain approximate redshifts by the flux density ratio between submm and radio wavebands 
(Carilli & Yun 1999, 2000), although uncertainties include an AGN contribution to the radio 
flux and the dust temperature in the galaxy, with a degeneracy between hotter galaxies at 
high redshift and cooler ones nearby (Blain 1999). 

• Relative importance of X-ray and UV stellar emission for reionization As previously 
noted, from a measurement of the radio synchrotron flux we can use equation (21) to estimate 
the level of inverse Compton X-ray emission which must be taking place. Redshifts and 
thus Ucmb(z) can be determined with Balmer line spectroscopy with NGST, while Ub can 
be estimated by assuming a magnetic field strength required to minimise the total energy 
density of the system (this is close to the value for energy equipartition between relativistic 
particles and magnetic fields) . Observations of local radio galaxies in non-thermal radio and 
X-ray emission yield field strengths consistent with the minimum energy value (Kaneda et 
al 1995). NGST can constrain the rest-frame UV emission, and a joint measurement of the 
rest-frame UV flux longward of Lya and the Balmer line flux (where Lh^ oc (1 — f e sc)^ion) 
can constrain f esc , the escape fraction of ionizing photons, where one roughly expects Jhu ~ 
40(1 — f esc ) Jir(R/ 1000). Thus, joint SKA/NGST observations could place limits on whether 
inverse Compton X-rays or stellar UV photons were energetically dominant and thus more 
important in reionizing the universe. 
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• Measuring magnetic fields in bright sources The very brightest sources out to z ~ 3 
will also be visible in X-ray emission with CXO; the fact that the same population of rela- 
tivistic electrons is responsible for non-thermal X-ray and radio emission can be confirmed by 
comparing respective spectral slopes. In this case, one can estimate the strength of magnetic 
fields, Ub ~ l ^ c Ucmb(z). Since dust obscuration is unimportant at both hard X-ray and 
radio wavelengths, the measurement should be fairly robust. 

• Distinguishing between synchrotron emission from AGN and supernovae To con- 
firm that such emission arises from star-forming regions rather than an AGN, one might look 
for signs of diffuse emission (note that the angular resolution of the SKA is ~ 0.1", while 
the angular scale of the virial radius of typical objects will be Q v i r ~ 

AGNs can be selected on the basis of color, as has been successfully carried out at lower 
redshifts (Fan 1999), using broad band NIR and MIR imaging with NGST. Finally, the line 
widths of the Ha, H/3 lines as observed with NGST may be considerably broader for an AGN 
(a ~ 1000 km s -1 ), due to line broadening by the accretion disc. 



4. Conclusions 

X-ray emission from a early star forming regions is predicted to be large and energetically 
comparable to UV emission. Non-thermal inverse Compton emission, which provides a good fit to 
local observations and should become increasingly important at high redshift, due to the evolution 
of the CMB energy density, is predicted to be the dominant source of X-rays. It introduces a whole 
host of physical consequences: the topology of reionization changes, becoming more homogeneous 
with much fuzzier delineation between ionized and neutral regions; reheating temperatures increase, 
with implications for feedback on structure formation and the observed width of Lya forest lines; 
the abundance of free electrons in dense regions increases, promoting gas phase H2 formation, 
cooling, and star formation. These effects will be considered in a companion paper (Oh 2000a). 
While direct detection of individual sources with CXO or Constellation-X appears difficult, we 
can hope to confirm the presence of relativistic electrons in high redshift objects by detecting 
non-thermal radio emission with the Square Kilometer Array. Given the CMB energy density at 
that epoch, this yields a minimal level of inverse Compton X-ray emission. Combined with NGST 
observations of rest frame UV emission, this will determine if stellar radiation or inverse Compton 
X-rays were the dominant factor in reheating and reionizing the universe. In addition, in this 
scenario a non-trivial fraction of the hard X-ray and gamma-ray background comes from inverse 
Compton emission at high redshift. In particular, it is possible to reproduce both the shape and 
amplitude of the gamma-ray background observed by EGRET, and predict that the majority of the 
smooth, isotropic gamma-ray background will remain unresolved by GLAST, which should display 
attenuation above ~ lOOGeV, from the pair production opacity due to ambient UV/IR radiation 
fields. 

Many of the conclusions in this paper depend upon a scenario in which the escape fraction of 
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UV ionizing photons is small, f esc < 10%. This assumption is well supported by observations in 
the local universe, as well as theoretical radiative transfer calculations of high redshift star forming 
regions, which predict very low escape fractions, f esc ~ 0.01 (Wood & Loeb 1999, Ricotti & Shull 
1999, although note that the latter authors predict a substantial escape fraction at low masses 
1O 7 M ). However, note that the latter ignore gas clumping and the multi-phase structure of the 
ISM. If for some reason the escape fraction is unexpectedly high (e.g. the supernovae blow holes 
in the ISM through which UV photons can escape), then the X-ray component is energetically 
subdominant and stellar UV radiation dominates the reionization of the universe. Even in this 
regime, the X-ray component still plays a role in heating the gas above 15 000 K, He II reionization, 
and promoting gas phase H2 formation, as these tasks cannot be accomplished by soft photons. 
The escape fraction of UV ionising photons in high redshift objects may eventually be deduced by 
comparing the IR and Ha fluxes observed by NGST (Oh 1999). 

The most uncertain aspect of this paper is the assumed level of X-ray luminosity, Lx ~ 
O.IEsn- I have calibrated the conversion rate via local X-ray and gamma-ray observations of 
starburst galaxies and individual supernova remnants within our galaxy, and argued that efficiency 
of all proposed X-ray production mechanisms either remains constant or increases with redshift. 
If the X-ray emission in local starbursts is primarily due to inverse Compton scattering of soft IR 
photons (Moran & Lehnart 1997, Moran, Lehnart & Helfand 1999), the empirical relation between 
star formation rate and X-ray luminosity (equation (1)) implies an electron acceleration efficiency 
of e ~ 10%, which lies at the upper limit of theoretical expectations. If instead the empirical 
relation (1) is approximately correct but the X-rays arise from a variety of emission mechanisms, 
the importance of X-rays for reionization still holds, but specific observational tests which rely on 
the inverse-Compton mechanism, such as the gamma-ray background observations (section (3.2)) 
and observations of radio synchrotron emission with the SKA (section (3.4)) will fail. 
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Fig. 1. — Comparison of stellar, QSO (with L v oc z^~ L8 ) and inverse Compton spectral energy 
distributions, all of which are normalised to have the same luminosity above 1 Ry. The inverse 
Compton SED is significantly harder than all others. At high densities (overdensities 5 ~ 10 4 
at z ~ 9), cooling of relativistic electrons by ionization losses becomes significant, causing the 
downturn of the spectrum at low energies. 
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Fig. 2. — Current best observational limits on the X-ray and gamma-ray background, against 
model predictions. The solid line represents model predictions for an IGM metallicity at z=3 of 
Z sa 1O _2 Z , while the upper and lower dotted lines represent generous limits of Z « 2.5 x 10 _2 Zq 
and Z ~ 1O _3 Z0 respectively. Note that inverse Compton radiation from high redshift can account 
for the majority of the observed gamma-ray background. 
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Fig. 3. — Number of objects with z > 5, z > 10 detectable in non-thermal emission in the 1° field of 
view of the SKA at 320 MHz; the limiting flux for a 10 a detection in a 10 5 s integration is displayed 

as a solid line. The dark solid line is for an efficiency factor f ra dio = f io^g ) (ot) ("oTT") = ^' 
the light solid line is for f ra dio = W~ 2 . Also shown as dotted lines is the detection rate and limiting 
flux at 4 GHz, when free-free emission is dominant, for a high efficiency star formation model 

(fstar = 17%). 



